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Abstract 

The radial motion along null geodesies in static charged black 
hole space-times, in particular, the Reissner-Nordstrom and stringy 
charged black holes are studied. We analyzed the properties of the 
effective potential. The circular photon orbits in these space-times 
are investigated. We found that the radius of circular photon orbits 
in both charged black holes are different and differ from that given 
in Schwarzschild space-time. We studied the physical effects of the 
gravitational field between two test particles in stringy charged black 
hole and compared the results with that given in Schwarzschild and 
Reissner-Nordstrom black holes. 
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1 Introduction 

The well known static, spherically symmetric black hole solutions in vacuum 
of Einstein's general relativity are given by the charged Reissner-Nordstrom 
and uncharged Schwarzschild solutions. In the non- vacuum case of Einstein's 
general relativity, several black hole solutions are known [1]- [3]. One of 
them is the stringy charged black hole discovered by Garfmkle, Horowitz and 
Strominger [3]. 

The study of timelike and null geodesies, the paths of freely moving par- 
ticles and photons, is the key to understand the physical importance of a 
given space-time. 

We wish to investigate in this paper the properties of the stringy charged 
black hole by studying its geodesic structure, that is, from the motion of 
photons. We compare the results by the aforementioned solutions in vacuum. 
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The space-time under consideration is almost identical to Schwarzschild 
space-time. The only differences are 

1. The gravitational energy, using M0ller's prescription [3], depends on 
the mass parameter M and on the charge Q (see Gad [5] ) , while in the 
Schwarzschild black hole is given only by the mass parameter M (see 
Xulu [6]). 

2. The areas of the spheres of constant r and t depend on the charge Q. 

The first aim of this paper is to sustain the second difference by analyzing 
the effective potential of radial motion along null geodesies. 

Compared to the Reissner-Nordstrom black hole, the stringy charged 
black hole exhibits several different properties [7J- [10]. For example, first, 
this solution has only one horizon at r = 2M and not two as is the case for 
the Reissner-Nordstronj^. Second, this solution has singularities at r = a, 
r = whereas the Reissner-Nordstrom has a singularity at r = 0. The ex- 
tremal solution occurs at Q 2 = 4M 2 e 2 * rather than at E 2 = M 2 as for 
Reissner-Nordstrom solution. 

It is therefore worthwhile to investigate other properties of the stringy charged 
black hole to see how this differs from the vacuum solutions. 

The equation of geodesic deviation gives the relative accelerations be- 
tween free test particles falling in a gravitational field and is a cornerstone 
to the understanding of the physical effects of the gravitational field and 
the geometry of space-time [9] . Geodesies deviations in vacuum space-times, 
namely, Schwarzschild and Reissner-Nordstrom space-times, are rigorously 
studied in [11]- [17]. Some studies on geodesic deviation can be found in [T8j - 
[20] . Ghosh and Kar [21] investigated geodesic motion and geodesic deviation 
in warped space-times with a a time dependent extra dimension. We do fol- 
low here the same approach in [16|,|T7J to obtain the relative accelerations of 
nearby test particles in stringy charged black hole. This is the second aim of 
this paper. 

The structure of the paper is as follows: In the next section, we investigate 
the geodesic curves and the circular photon orbits in the stringy charged black 
hole. In section 3, we study the photon trajectories in Reissner-Nordstorm 
by analyzing the properties of effective potential. Section 4 concerns with 
the tidal forces between nearby test particles in stringy charged black hole. 
In section 5, we investigate the behavior of geodesic deviation vector along 
a timelike geodesic near the singularities r = and r = a, by considering a 

2 The Reissner-Nordstrom has two horizons given by the quadratic equation r 2 — 2Mr + 
e 2 = 0. 
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special case. Finally, in section 6, brief summary of results and concluding 
remarks are presented. 

Through out this paper Latin indices run from to 3; Greek indices run from 
1 to 3. Units are such that G = 1 = c, (G is the gravitational constant; c is 
the velocity of light). 



2 The equations of geodesic motion 

In recent years there is considerable interest in obtaining black hole solu- 
tions in string theory and investigating their properties (see references [3], [7] 
and [22]). There are two metrics in this theory, which are called the sigma- 
model (or string) metric and Einstein metric. For uncharged static, spheri- 
cally symmetric black hole, the solution in the low energy is the same as the 
Schwarzschild solution. This is, however, not the case when the black hole is 
charged. 

We are interested to investigate the geodesic curves of a static spherically 
symmetric charged black hole. The line element representing this space-time 
is given by Garfmkle et al [3]. 

ds 2 = _ (1 _ ™ )dt 2 + (1 _ ^yl dr 2 + (1 _ ^) r 2( d 02 + sin 2 #2)j (2 -Q 

ty ry ry 

where 

a = ^-exp(-2%) 1 (2.2) 

M and Q are, respectively, mass and charge parameters; $o is the asymptotic 
value of dilaton field. 

The equations and constraint for geodesies are given as 

x a + T a bc x b x c = 0, 

where a superposed dot stands for a derivative with respect to the affine 
parameter r associated to the geodesic, x a are the coordinates of a space- 
time point on the geodesic and e = — 1 or 0, for timelike or null geodesies, 
respectively. 

The geodesic equations for the line element (12. ip are given by 

2M , s 

t+— — rt = 0, 2.3 

r(r-2M) ' v ' 



3 



M . 2 (r - 2M){a - 2r) A2 r - 2M a - 2r . 2 - 2 

r ; -r 2 + - - -9 2 + - - sin 2 6»0 2 + 

r(r-2M) 2r 2r Y 

M{r ~ 2M) P = 0, (2.4) 

+ ^ - 2r j rfl - sin 9 cos 6>0 2 = 0, (2.5) 
r(a — r) 

(a — 2r) ■ 

4> + ^ fr</> + 2cot000 = O, (2.6) 

r(a — r) 

The constraint of timelike or null geodesies for the line element (12. ip is 
given by 

_ (1 _ ™)i 2 + (1 - — )- l f 2 + (1 - -)r 2 {9 2 + sin 2 6tf) = e. (2.7) 

ry> ty f 

Let a geodesic 7 be given by j(t) = (t(r),r(r), 0(r), 0(r)). We suppose 
without loss of generality that 7(r ) = (t(r ), r(r ), |, 0(ro)) and for all r : 
= f (equatorial orbits ) that is the particle has at start, and continues 
to have, zero momentum in the ^-direction; thus 9 = 0. Consequently, the 
equation of geodesies (I2.3P and (12. 6p have the straightforward first integrals 

i = *((^fij). (2 ' 8) 

i = -r^— v (2.9) 
r[r — a) 

The integration constants c\ and C2 can be found, if we know the initial 
conditions 7(r ) and ^-(to) for some r G 3?. In the case of Schwarzschild and 
Reissner-Nordstrom space-times, Clarke [23] and also Wald [21] demonstrate 
that Ci represents the total energy, E, per unit rest mass of a particle as 
measured by a static observer, and C2 represents the angular momentum, 
L, per unit mass of a particle (see also [LT]). We recognize the constant 
s to represent the rest energy per unit mass for massive particles (timelike 
curve, e = — 1) or the rest energy for massless particles (null curves, e = 0), 
traveling along the given geodesic [TT][T2] . 

Using equations 02 .8p . (12.91) and 9 = | in the constraint equation (12. 7p . 
we have 

r 2 = £ 2 -(l-— )[-^ + T i -^]' ( 2 ' 1Q ) 
r L r(r — a) 

Notice that we have not used equation (12 .4p for the following reason: If we 
substitute r and its derivative with respect to r from the constraint equation 
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(I2.7P and using equations (12 .8p and f 12 . 9 j) in equation (12. 4p . then it is satisfied 
identically. 

Equation (I2.10p can be written as 

r 2 = E 2 — V 2 (r), (2.11) 

where V(r) is the "effective potential" defined by 

V 2 (r) = (l-— )[-^+ T i -^]' ^ 
r L r(r — a) J 

Equation (12. lip is in the form of the equation of a one-dimensional problem 
for a particle in a potential field V(r). 

Since the left side of equation (12.101) is positive or zero, the energy E of 
the trajectory must not be less than the potential V. So for an orbit of a 
given E, the radial range is restricted to those radii for which V is smaller 
than E. 

In the following we investigate the circular photon orbits in the stringy 
charged black hole by analyzing the properties of the effective potential. 
In the case of the photon trajectories, putting e — 0, the effective potential 
(gUZD takes the form 

^ 2 (r) = (l-— )[ T ^]- (2.13) 
r L r(r — a) J 

Differentiating equation (12.111) with respect to r gives 

2(^ r )(^ 2r ) dV 2 {r) dr 
^dr dr 2 dr dr' 

or 

It is clear from equation (I2.14p that a circular orbit (r= constant) is possible 
only at a minimum or maximum of V 2 (r). We can quantitatively evaluate 



dr 



:i-^)[-r^] 

r r(r — ct) 



and get 

r = - ((6M + a) ± ^(eM + a) 2 - 32Ma) . (2.15) 

Since a must be less than r, to keep the line element (12.11) to be in Lorentzian 
metric, then for a = 2M and a > 2M no circular is possible. From the 
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relation (12. 2p a is smaller than M (a < M), therefor the larger of two roots 
given by equation (12.151) locates the maximum of the potential-energy curve 
V 2 (r) defined by equation (12.131) . Consequently, the unstable circular orbit 



is always at the same radius r = \ [{6M + a) + y/(6M + a) 2 - 32Maj , 

regardless of L. For the smaller root no circular is possible. 

There are many possible ways to compare the effective potential for a 
massless particle in stringy charged black hole with that in Schwarzschild 
black hole. Here we assume a < M and fix the value of angular momentum 
L; the values of r and a are calculated to obtain the plot of the effective 
potential. The choice of a << M allows the two paths to be the same. (This 
comparison is shown in figures (1) - (5).) 

We notice that, when a = 0, that is the Schwarzschild case, the radius is 
r = 3M which is the same radius obtained by Schutz [13] in the Schwarzschild 
black hole. 

3 Reissner-Nordstrom metric 

A well-known simple solution of the Einstein-Maxwell equation is the Reissner- 
Nordstrom solution. This solution represents a non-rotating charged black 
hole. Discussions of the basic properties of this solution can be found in many 
places including works by Chandrasekhar [T3], and Hawking and Ellis [IB] . 
The metric is defined on a four- dimensional manifold and its typically written 
in the form 



ds 2 = _(1 _ ™ + U dt 2 + (1 _ *m + U-l dr 2 + r 2 {M 2 + gin 2 #2) (31) 



where m represents the gravitational mass and e the electric charge of the 
body. 

The equation of geodesies and the constraint of geodesic for the line ele- 
ment (13.11) are given as 




2 



2 



t + 



2(-m + f) 



ft = 0, 



(3.2) 



r 



2 — 2mr + e 2 




r 



4>+ -rj> + 2 cot 00<j> = 0, (3.4) 



r 




(1 - — + ^)" 1 f 2 + r\Q 2 + sin 2 9tf) = e. (3.5) 
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We will assume, as section 2, that the orbit is in the 9 = | plane. Equation 
(13.31) shows that if 9 — | and 9 = initially, then 9 = and the orbit remains 
in this plane. 

Equations (I3.2p and (13 .4[) can be integrated directly, giving 



r 2 — 2mr + e 2 

; c 4 



(3.6) 



(3.7) 



where the integrating constant C3 represents the energy, E, (at r — > 00) of a 
test particle and C4 the angular momentum, L. 

Substituting (13. 6 p and (13. 7p in equation (13 .5p and using the condition 
9 = 0, we get 

2m e 2 . r L 2 



r 2 = E z -(l- 



Equation (13. 8 p can be written as 



+ 



[- £ + -f] 



r 2 = E 2 -V 2 (r) 



where V(r) is the "effective potential" defined by 

-,2 



V 2 (r) 



2m e\r L 2 , 



(3. 



(3.9) 



(3.10) 



In this paper we restrict our tension to the photon trajectories, by putting 
e = in equation (13.101) . the effective potential becomes 



V 2 (r) 



2m e 2 L 2 



[».2 J " 



Differentiating (I3.9P with respect to r, as in the previous section, we get 

d 2 r 1 d 



(3.11) 



— = -- — (V 2 (r)). 
dr 2 2dr K 1 JJ 



(3.12) 



This equation shows that a circular orbit (r= const.) is possible only at a 
minimum or maximum of V 2 (r). We can quantitative by evaluating 



d 

dr 



^ 2m e 2 r L 2 



which gives 



3m 



1 ± a/1 - 



9m 2 



(3.13) 



This equation shows that the two radii are identical for — = (in this case 
a minimum radius for photon tmin = ^p) and do not exist at all for — > ^p. 
This indicates a qualitative change in the shape of the curve of V 2 for small 
— . For — < the larger of the two roots given by equation (I3.13P locates 
the minimum of the potential- energy curve V 2 defined by equation (13.111) . 
while the smaller root locates the maximum of the potential-energy curve. 
Therefore, the circular orbit of the larger radius will be stable in contrast to 
circular orbit of the smaller radius which will be unstable. 



4 Geodesic Deviation 

In this section, we use the tidal forces between free test particles falling in 
a gravitational field to investigate different properties between the stringy 
charged and vacuum black holes. 

Consider a sphere of two non-interacting particles falling freely towards 
the center of the Earth. Each particle moves on a straight line, but nearer 
the Earth fall faster because the gravitational attraction is stronger. This 
means that the sphere does not remain a sphere but is distorted into an 
ellipsoid with the same volume. The same effect occurs in a body falling 
towards a spherical object in general relativity, but if the object is a black 
hole the effect becomes infinite as the singularity is reached. Jacobi vector 
fields provide the connection between the behavior of nearby particles and 
curvature, via the equation of geodesic deviation (Jacobi equation) 

D 2 n a 

u 'I to „.b„.c„d 



Dt 2 



+ Rf cd v b vY = 0, (4.1) 



where v a are the components of the tangent vector to geodesic and if are 
the components of the connecting vector between two neighboring geodesies. 

In order to investigate in detail the behavior of Jacobi fields we consider a 
congruence of timelike geodesies (path of particles) with timelike unit tangent 
vector v (g(v,v) = —1). We define at some point q on the geodesic 7(7-) 
dual bases eg, e^, e^, eg and e°, e^, e^, e„ of the tangent space T q M and dual 
tangent T*M respectively in the following way [16] : We choose eg to be v a 
and e", e%, eg as unit spacelike vectors, orthogonal to each other and to v a . 
If we parallelly propagate the basis along the timelike geodesic 7(7") (that is, 
Tj^eg, = 0, a = 1, 2, 3), eg will remain equal v a , and e", eg, eg will remain to 
orthogonal to v a (see [15] p. 80). The frame eg, e^, eg, eg is called "parallel 
transported" (PT) frame. The orthogonal connecting vector, r/ a , between two 
neighboring timelike geodesies may be expressed as rf = rj a e^, (rf = e^,r] a = 
0). 



s 



The geodesic deviation vector rf satisfy the following equation 

Drf 



Dt 



rf a v\ (4.2) 



2„a 



D 2 7] 



+ R^> h v c e*rf = 0, (4.3) 

where rf are the space-like components of the orthogonal connecting vector 
rf connecting two neighboring particles in free fall; rf = 0. The tilde denotes 
components in the PT frame and the components of the Riemann tensor R% dc 
are given by 

Rlc = <e{ele h d W fgh . (4.4) 
From (12. ip the frame e£ in Reissner-Nordstrom metric is given by : 



e a =(l-^)-*(0, 0,0,1), 
e a = (i_M)i(i j0j 0,0), 

e a 2 = {J -?—(0, 1,0,0), 1 J 



3 =^(0-0,1,0). 



The components of rf can be written as follows 

Using (I4.4p . (14. 5p . v a = eg and the components of Riemann tensor for the 
metric (12. ip (see appendix), in (14. 3p . we get 

£>V _ 2Af„r 

Dr 2 ~~ W'l ' 

DV _ M(a-2r) fl , , fi s 

£>V _ M(a-2r) A, 
Dt 2 2r 3 (r-a) " ' 



In order to write equation (14. 6 p in terms of ordinary derivative, we must 

Dr 2 



evaluate the second covariant derivative -J^. Using eg = v a , equation (14. 2 p 



takes the form 

D,f A '° + f^V, (4.7) 



Dt dr 
where = e™e f a e 9 b F e fg . 

Differentiating (14. 7p covariantly and using the Christoffel 
components of metric (12 .ip (see appendix), we can write 
(14.61) in the form 
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if, (4.8) 



dr 2 r 3 
<Prf M(a-2r) e 



dr 2 2r 3 (r — a) 
d 2 rf M(a - 2r) 



^, (4.9) 
rf. (4.10) 



<ir 2 2r 3 (r — a) 

Equation (14.81) indicates that tidal force in radial direction will stretch 
an observer falling in this fluid. To keep the line element (12. ip to be in 
Lorentzian metric, a should be less than r. Therefore equations (14. 9p and 
(I4.10p indicate a pressure or compression in the transverse directions. 



5 Solution of the equations (14.81) - ( 14.101) 



To solve the equations (I4.8p - (l4.10p . we consider the special case when freely 
falling particles have zero angular momentum {L = 0). From equation (I2.10p 
we obtain the relation between the radial coordinate, r, and the affine pa- 
rameter, r, in the following form 



dr 
d^ 



(E 2 



D + 



2M 



(5.11) 



Using this relation in equations (I4.8p - (l4.10p . we consider, without loss of 
generality, the case when E 2 = 1, we obtain 



d 2 rf Idrf _ 1 
dr 2 r dr r 

d 2 rf 1 drf 
dr 2 2r dr 

d 2 rf 1 drf 



(2r - a) e 

'/ ; 



dr" 2 



r dr 



Ar 2 (r — a) 
4r 2 (r — a) 



(5.12) 

(5.13) 
(5.14) 



Solving the above equations, using MAPLE, we obtain the space-like com- 
ponents of the geodesic deviation vector in the following form 



+ C 2 r" 



(5.15) 



C 3 F 1 (l + b,b,-l- + 2b, -)r b (r-a) + 
z z a 
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C A F 2 {1 - 6,2 - 6, | ~ 26, -)r^- b \r - a), (5.16) 
Z I a 

rf = C.F, (i + 6, 6, -i + 26, -)r b (r - a) + 

C 6 F 2 f| - 6,2 - 6, | - 26, ^)r(f- fc )(r - a), (5.17) 

where Ci, C2, C3 and C4 are constants of integration and 6 = | ± Fi and 
F 2 are hypergeometric functions. 

The spacial solution given by (I5.15p -f l5.17l) allows one, by a suitable choice of 
initial conditions and constants appearing in f l4.8p - fl4.10p . to describe various 
physical situations connected with relative motion of freely falling particles. 

At the singularity r = the radial component rf of geodesic deviation 
vector becomes infinite, while the transverse components rf and rf vanish. 
Therefore, the behavior of geodesic deviation vector in the space-time under 
consideration is the same as in the case of Schwarzschild space-time, when the 
singularity, r = 0, is approached. In fact the Schwarzschild singularity has 
time-like geodesic for which one component of the geodesic deviation vector 
becomes infinite while the other two vanish [11]. The singularity r = 
is a strong curvature singularity, as the Schwarzschild singularity [26], for 
the two vanishing components of geodesic deviation vector go to zero faster 
than the remaining component becomes infinite. The misbehavior of these 
components could in concert so that it is not visible in the volume element 
magnitude. 

At the singularity r = a, the component rf has a finite value, while the 
other two components, rf and rf ', vanish. Then the volume element, defined 
by the three space-like components of geodesic deviation vector [25], [26] . 
along a time-like geodesic vanishes as the geodesic approaches this singularity. 

6 Conclusion 

In this paper we have studied the circular photon orbits in charged black holes 
by analyzing the properties of effective potential. Considering the light-like 
geodesies, we classified and analyzed the different cases between the stringy 
charged black hole and the vacuum solutions. These differences arose from 
considering the orbits associated with stable and unstable circular orbits. In 
the context of the Schwarzschild geometry there is an unstable circular orbit 
which is always at the same radius, r = 3M. In the Reissner-Nordstrom 

3 A singular point is called a strong curvature singularity if any object hitting it is 
crushed to zero volume |26j . 
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space-time there are two radii, the circular orbit of the larger radius will be 
stable while that of the smaller radius will be unstable. In the case of stringy 
charged black hole there is an unstable circular orbit as the Schwarzschild 
black hole, but the difference is that the radius in the case of stringy charged 
black hole depends on the charge Q. 

Equations fl4.8p - fl4.10p provide the explicit expressions of the relative ac- 
celerations in a stringy charged black hole. Two comments are worth making 
about expression f!4.8p . First there is no divergence in the radial direction at 
r = 2m. Secondly, the tidal field at the horizon in radial direction is larger 
for smaller black hole. This is simply because 

d 2 V r 2M r 1 . x 

77 atr ~ M - 

The radial component of the geodesic deviation vector field, equation 
(14. 8p . is the same as the component obtained in the case of Schwarzschild 
black hole and different from the component obtained in the case of Reissner- 
Norstrom. Consequently, in the case of stringy charged and Schwarzschild 
black holes, the tidal forces in the radial direction will stretch an observer 
falling in these black holes, while in the case of Reissner-Norstrom the ra- 
dial component depend on the quantity 2e 2 — 2mr to indicate a tension or 
stretching in the radial direction. In the transverse directions the relative ac- 
celeration components (see equations (14 .9p and (I4.10p ) indicate a compression 
in these directions. This property is similar as in the case of Schwarzschild 
black hole, while in the case of Reissner-Nordstrom depends on the quantity 
e 2 — 2mr to indicate compression or tension in these directions. 

To describe various physical situations connected with relative motion of 
freely falling photons, in stringy charged black hole, we tried to solve equa- 
tions (I4.8p - fl4.10p by considering a special case. We found that any object 
hitting the singularity r = a is crushed to zero volume. 



Appendix 

We U.SG (x ■ X , X . X 3 ) = (t,r,9,(p) so that the non- vanishing Christoffel sym- 
bols of the second kind of the line element fl2.ll)') are 



-pl M p2 a-2r 

1 11 — r(2M-r)' 1 12 ~~ 2r(a-r) ' 

r l _ {r-2M)(a-2r p 3 _ q-2r 



22 2r ' 1 13 — 2r(a-r) ' 

= (r-2A0 (a -2r) ^2 g> ^ = _ sin ^ cos ^ 

Too r| 3 =cottf, 

■pO M 

1 10 r(r-2M) ' 
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The non-zero Riemann tensor are: 



n 212 — 4r 2 (a-r) ' 

pl _ 2M(r-a)(2r-a)+a 2 (2M-r) - 9 a 

n 313 - 4r 2 (a-r) S111 P ' 

pi _ 2M(2A/-r) 

-n-OlO ~~ r 4 j 

p2 _ 8Mr(r-a)-a 2 (r-2M) ^2 g 

-"-323 ~~ 4r 2 (r-a) b ' 

D2 _ R3 _ A/(o:-2r)(2A/-r) 

-"-020 — -"-030 — 2r 4 (r-a) 

References 

J. D. Bekenstein, Ann. Phys. (N. Y.), 91, 75 (1975). 

E. Ayon-Beato and A. Garcia, Gen. Relat. Grav. 31, 629 (1999). 

D. Garfinkle, G. T. Horowitz and A. Strominger, Phys. Rev. D43, 3140 
(1991); D45, 3888(E) (1992). 

C. M0ller, Ann. Phys. (N. Y.), 4, 347 (1958). 

R. M. Gad, Astrophys. Space Sci. 295, 459 (2004). 

S. S. Xulu, Astrophys. Space Sci. 283, 23 (2003). 

J. Preskill, P. Schwarz, A. Shapere, S. Trivedi and F. Wilczek, Mod. 
Phys. Lett., 6, 2353 (1991). 

J. A. Harvey and A. Strominger, "Quantum aspects of black holes", 



Preprint EFI-92-41, [hep-th/ 9209055 



F. A. E. Piran, Phys. Rev. 105, 1089 (1957). 

K. P. Tod, Proc. R. Soc. Lond. A388, 467 (1983). 

C. Misner, K. Thorne and J. Wheeler, (1973), "Gravitation", Freeman, 
San Francisco. 

R. Adler, M. Bazin and M. Schiffer, (1975), "Introduction to General 
Relativity", (McGraw-Hill, New York, 2nd. ed.). 

B. F. Schutz, (1985), " A First Course in General Relativity" ( Cambridge 
Uni. Press, Cambridge, London, New York, New Rochelle, Melbourne 
Sydney). 



13 



[14 

[15 

[16 

[17 

[18 

[19 

[20 
[21 

[22 

[23 

[24 
[25 



S. Chandrasekhar, (1983), "The Mathematical Theory of Black Holes", 
(Oxford Uni. Press, Cambridge, England). 

S. W. Hawking and G. F. R. Ellis, "The Larger Scale Structure of Space- 
time", (Cambridge Uni. Press, Cambridige). 

R. D'Invermo, "Introducing Einstein's Relativity", Oxford University 
Press, New York, (1992). 

M. Abdel-Megied and R. M. Gad, Chaos, Solitons and Fractals, 23, 313 
(2005). 

R. Kerner, J. Martin, S. Mignemi and J. E. Van Holten, Phys. Rev. D 
63, 027502 (2001). 

G. F. R. Ellis and H. van Elst, "Deviation of geodesies in FLRW space- 



time geometries", arXiv: gr-qc/9709060 



R. Koley, S. Pal and S. Kar, Am. j. Phys., 71, 1037 (2003) 

S. Ghosh and Sayan Kar, "Geodesies and geodesic deviation 
in warped spacetime with a time dependent extra dimension", 
larXiv:0904.2321^ 1 [gr-qc] . 

A. Shapere, S. Trivedi and F. Wilczek, Mod. Phys. Lett. A6, 2677 
(1991). 

C. J. S. Clark, (1979), "Elementary General relativity", (Edwaed Arnold, 
London) . 

R. M. Wald, (1984), "General Relativity", (Chicago and London). 

R. Penrose, (1972), "Techniques of Differential Topology in Relativity", 
(SIAM, Philadelphia, p. 60). 



[26] F. J. Tipler, Phys. Lett. 64A, 8 (1977). 



11 



